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Abstract 

The classical definition of Prym varieties deals with the unramified 
covers of curves. The aim of the present paper is to give explicit alge- 
braic descriptions of the Prym varieties associated to ramified double 
covers of algebraic curves. We make a careful study of the connection 
with the concept of algebraic completely integrable systems and we 
apply the methods to some problems of Mathematical Physics. 
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1 Introduction 

During the last decades, algebraic geometry has become a tool for solving 
differential equations and spectral questions of mechanics and mathematical 
physics. The present paper consists of two separate but related topics : the 
first part purely algebraic-geometric, the second one on Prym varieties in 
algebraic integrability. 

Prym variety Pr ym(C/Co) is a subabelian variety of the jacobian variety 

Jac(C) = Pic°(C), 

= H l (O c )/H l (C,Z), 

constructed from a double cover C of a curve Cq: if a is the involution on C 
interchanging sheets, then a extends by linearity to a map 

a : Jac(C) — > Jac(C). 

Up to some points of order two, Jac(C) splits into an even part and an odd 
part : the even part is Jac(Co) and the odd part is a Prym(C/Co). The clas- 
sical definition of Prym varieties deals with the unramified double covering 
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of curves and was introduced by W. Schottky and H. W. E. Jung in relation 
with the Schottky problem [6] of characterizing jacobian varieties among all 
principally polarized abelian varieties (an abelian variety is a complex torus 
that can be embedded into projective space). The theory of Prym varieties 
was dormant for a long time, until revived by D. Mumford around 1970. It 
now plays a substantial role in some contemporary theories, for example in- 
tegrable systems [1,2,5,8,10,14,15,17], the Kadomtsev-Petviashvili equation( 
KP equation), in the deformation theory of two-dimensional Schrodinger op- 
erators [19], in relation to Calabi-Yau three-folds and string theory, in the 
study of the generalized theta divisors on the moduli spaces of stable vector 
bundles over an algebraic curve [7,11],... 

Integrable hamiltonian systems are nonlinear ordinary differential equa- 
tions described by a hamiltonian function and possessing sufficiently many 
independent constants of motion in involution. By the Arnold-Liouville the- 
orem [4,16], the regular compact level manifolds defined by the intersection 
of the constants of motion are diffeomorphic to a real torus on which the 
motion is quasi-periodic as a consequence of the following differential geo- 
metric fact; a compact and connected n-dimensional manifold on which there 
exist n vector fields which commute and are independent at every point is 
diffeomorphic to an n-dimensional real torus and each vector field will define 
a linear flow there. A dynamical system is algebraic completely integrable 
(in the sense of Adler-van Moerbeke [1]) if it can be linearized on a complex 
algebraic torus C n /lattice (=abelian variety). The invariants (often called 
first integrals or constants) of the motion are polynomials and the phase 
space coordinates (or some algebraic functions of these) restricted to a com- 
plex invariant variety defined by putting these invariants equals to generic 
constants, are meromorphic functions on an abelian variety. Moreover, in 
the coordinates of this abelian variety, the flows (run with complex time) 
generated by the constants of the motion are straight lines. 

One of the remarkable developments in recent mathematics is the inter- 
play between algebraic completely integrable systems and Prym varieties. 
The period of these Prym varieties provide the exact periods of the motion 
in terms of explicit abelian integrals. The aim of the first part of the present 
paper is to give explicit algebraic descriptions of the Prym varieties associ- 
ated to ramified double covers of algebraic curves. The basic algebraic tools 
are known and can be found in the book by Arbarello, Cornalba, Griffiths, 
Harris [3] and in Mumford's paper [18]. In the second part of the paper, we 
make a careful study of the connection with the concept of algebraic com- 
pletely integrable systems and we apply the methods to some problems such 
as the Henon-Heiles system, the Kowalewski rigid body motion and Kirch- 
hoff 's equations of motion of a solid in an ideal fluid. The motivation was the 
excellent Haine's paper [10] on the integration of the Euler- Arnold equations 
associated to a class of geodesic flow on SO(A) . The Kowalewski's top and 
the Clebsch's case of Kirchhoff 's equations describing the motion of a solid 
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body in a perfect fluid were integrated in terms of genus two hyperelliptic 
functions by Kowalewski [13] and Kotter [12] as a result of complicated and 
mysterious computations. The concept of algebraic complete integrability 
(Adler-van Moerbeke) throws a completely new light on these two systems. 
Namely, in both cases (see [14] for Kowalewski's top and [10] for geodesic 
flow on 50(4)), the affine varieties in C 6 obtained by intersecting the four 
polynomial invariants of the flow are affine parts of Prym varieties of genus 3 
curves which are double cover of elliptic curves. Such Prym varieties are not 
principally polarized and so they are not isomorphic but only isogenous to 
Jacobi varieties of genus two hyperelliptic curves. This was a total surprise 
as it was generically believed that only jacobians would appear as invariants 
manifolds of such systems. The method that is used for reveiling the Pryms 
is due to L.Haine [10]. By now many algebraic completely integrable systems 
are known to linearize on Prym varieties. 

2 Prym varieties 

Let 

ip : C — » C , 

be a double covering with 2n branch points where C and Cq are nonsingular 
complete curves. Let 

cr : C — >C, 

be the involution exchanging sheets of C over Co = C/a. By Hurwitz's for- 
mula, the genus of C is g = 2go + n — 1, where go is the genus of Co- Let 
Oc be the sheaf of holomorphic functions on C. Let S 9 (C) be the g-th sym- 
metric power of C (the totality of unordered sets of points of C). On S 9 (C), 
two divisors T>\ and T>i are linearly equivalent (in short, T>\ =T>2) if their 
difference T>\ — T>2 is the divisor of a meromorphic function or equivalently 
if and only if 

where is the sheaf of holomorphic 1-forms on C and 7 is a closed path on 
C. We define the jacobian (Jacobi variety) of C, to be 

Jac(C) =S 9 (C)/ = . 

To be precise, the jacobian Jac(C) = Pic°(C) of C is the connected compo- 
nent of its Picard group parametrizing degree invertible sheaves. It is a 
compact commutative algebraic group, i.e., a complex torus. Indeed, from 
the fundamental exponential sequence, we get an isomorphism 

Jac(C) = Pic°(C), 

~ fl°(c,^)7fri(c,z), 

~ C 9 /Z 29 , 
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via the duality given by Abel's theorem. Consider the following mapping 
S 9 ( C ) c°/Ln, J2^^J2J (v 1 ,...,u g )=t(k 1 ,...,kg), 



fc=l 



k=l 



where (cji, . . . ,uj g ) is a basis of f^, Lf2 is the lattice associated to the period 
matrix Q and fi±, . . . ,fi g some appropriate variables defined on a non empty 
Zariski open set. Let (ai, ...,a go , ...,a g ,bi, ...,b go , ...,b g ), be a canonical ho- 
mology basis of i?i(C,Z) such that 

a (en) = a fl0+n , cr (a fl0 ) = a 9 , 

<7 (%)+l) = -Ogo+1, —,Cr {a go +n-l) = —Ugo+n-l, 

° ( fe i) = b go+n, -, v {bg ) = bg, 

° (bg +i) = -bg +i, ...,a(b go+n -i) = -b go+n -i, 

for the involution a. Notice that ip (a go+1 ) , ip (a 9o+n _i) , p (b go+1 ) , p (6 go+n _i) 
are homologous to zero on Co- Let . . . ,co g ) be a basis of holomorphic 
differentials on C where to go+n , . . . ,u g , are holomorphic differentials on Cq 
and 

-uj, l<j<go + n-l, 
Vj, go + n < j < g, 

the pullback of Uj. The period matrix $7 of Jac(C) is explicitly given by 



<T*{Uj) 



Q = 



A B C D E F 
G H I J K L 



where A, . . . ,L denote the following matrices : 



A 



y iai ^90+™-! 



B 



c 



( In 

Ja go+ i 



\ Ia gQ+1 W S0+n-l 
/ L + ^ 

Ja 90+ n 



Ja Sn+n-l 



(1) 



I? 



bi 



lb w l \ 



Jb, 



bgo ^90+n-l / 
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and 



E = 



6 S0 + 1 



'so+ r 



\ 



V ^> 90 + l a; 90+"- 1 ' ' ' -^flo+n-l W SO+n-l / 
/ /ft ^ ^1 • • • /ft ^1 \ 



G 



# = 



J 



X = 



L = 



Notice that 



U>4 



\ /6 fl o+« • • • Ib g ^90+n-l J 

I Li W SO+n • • • J a ^g +n \ 



\ /ai W J ••• /a 9(J ^9 y 
/ /a 90+1 W 90+n • • • /a so+n _i U 90+n \ 



la 



Ia a U 90+n \ 



) 



/ /ftj u 9Q+n ■ ■ ■ k V go+ n \ 



\ fbi U 9 ■■■ k go U 9 ) 
//ft 



/ft 90+ „-i ^9 J 
k a U 90+n \ 



a 90 + l 



T ( a 9() + l) 



/ 

J a an + l 



'J: 



1 < j <9o + n-l, 



Ja gQ+ i 

-/a 90+1 ^' 5o + n<J<5, 



5 



/ 

I > 



Uj 



L uj, l<j< go + n-l, 
L ^ , w i> 9o + n<j <g, 

90+n-l 



CT ( 6 SO + l) 

a* {uj) , 



90 + 1 



/ 6fl0+1 w i> l<j<9o + n-l, 



-Jb go+1 u i' 9o + n<j<g, 



' fe 9 + l 



/ 

Jb g . 

J a„ 



Uj 



Uj 



a 9Q + " 



4 o+1 ^-, l<i<5o + n-l, 
4 0+ i^' 9o + n<j<g, 



I 



Uj, 



cr(ai) 
J ai 



and 



/ 



Uj 



a 90+ n 



-/„ Wj, l<i< go + n-l, 



"90 



L 90 w i> go + n<j<g, 



/ = 



U 



.1 ' 



<r(6i) 



/^Wj, l<J<5o + n-l, 
/ 6l Wj, go + n<j<g, 



/ 

Jb„ 



Uj = 



j bgo UJj, l<j< go + n-l, 
Ib go u ii 9o + n<j<g. 
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Then, 

C= -A, F = -D, H = 0, I = G, K = 0, L = J, 

where O is the null matrix and therefore 

f A B -A D E -D\ 
\G G JO J J' 

= ( C\ C2 C3 C4 C5 C% ) . 
By elementary column operation, we obtain the following matrices : 

fii = ( C\ C2 C\ + C3 C4 C5 C4 + Cq ) , 

f A B O D E O \ 
\G O 2G J O 2J J ' 



Q2 — ( Ci C2 Ci — C3 C4 C5 C4 — 6*6 ) , 

/ 4 5 24 D E 2D \ 
\G O O J O O ) ' 



$^3 — ( C\ — C3 C2 C4 — Cq C5 Ci + C3 C4 + Cq ) 

2A B 2D E O O 
O O O O 2G 2J 

r O 
O 2A 



where 
A 



(G J), 



V 



Ia go J bl ■ ■ ■ J bgQ ^g +n 



Ja 



9(1 



k ^9 



b ^9 
"so y 
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and 



r = (2A B 2D E ) , 



(3) 



2 Jai W 90+n-l \ T 



2 /a W l 

Ja 30 
Ja 9 + 1 



2 Ja 90 Vgo+n-1 
Ia go+1 



• • io so+ „_i W flo+n-l 
2 Jfcx w flo+n-l 



2 Ib gQ W l • • • 2 J& on ^90+n-l 



90 



Let 



in 



9 , 
i=l Ja * 



I "1 \ 



ui 



: mi,m £ Z}, 



be the period lattice associated to f2. Let us denote also by Ln 15 £f} 2 , £n 3 
and La the period lattices associated respectively to fii, f^, and A. Since 



— Lq 1 — Ln 2 — Lq 3 , 



it follows that the maps 



C 9»+«-l/ Lr 



\ tgo+n-l / 



modLr ^ C 9 /L a : 



( ti \ 



V J 
( \ 



modLn, 



C 9 " /2L A : 



mod2L A ^ C 9 /L n : 





ti 

\ ho J 



modL^, 
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are injectives. Therefore, the tori C 9o+n 1 / 'Lp and C 9o /2L 1 \ can be embed- 
ded into C 9 1 Lq and the map 

C»/Ln 3 = C» 0+n - 1 /^r e C3°/2L A — > C^/L^, 



modL^., 



V h ) 



modLn, 



shows that the jacobian variety Jac (Co) intersects the torus C 9o+n 1 / Lr in 
2 2 9o points. We have then the following diagramm 



C 

c 

if 

Co 



u 



KerN v 
I 

Jac(C) 

1% 

Jac(Co) 

if* 
Jac(C) 

i<y 

Jac(C) 



— Jac*(C) 

in; 

^- Jac* (C ) 
C H\C,C%) 

Jac(Co) 



N, 



H\Co,{<p*OcY) 
H\Co,0* Cq ) 



where 



u : z i — ► divisor class(z —p), p € C, fixed, 
uo ■ zo i — > divisor class(zo — po), Po € Co, fixed withpo = <p(p), 
N v : Jac(C) — ► Jac(C ), ^rriiPi i — > ^mupjpj), 
is the norm mapping and 

Jac*(C) = dual of Jac(C). 

The norm map N v is surjective. Moreover, the dual Jac* (C) of Jac(C) is 
isomorphic to Jac(C). The Prym variety denoted Prym(C /Co) is defined by 

Prym{C/Co) = (h° (C, n c ) ") * /H, (C, Z)" , 
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where —denote the -1 eigenspace for a vector space on which j acts. Let 
D € Prym(C/C ), i.e., 



\ 


2 Jo U 90+n 

V 2/ % J 



rD 

Jo ^go+n 



G La, 



f 
Jo 



CJ 



rD 

Jo 

rD 

Jo W 90+« 



V /c 



W 9 



/ *1 \ 







/ V / 



Consequently, the Abel-Jacobi map 



( lo 



Jac(C) -^C 9 /L n , D 



rD 

Jo w 90 + «-l 
rD 

Jo ^go+n 
\ Jo u g ) 

maps Prym{C / Cq)) bi-analtically onto C 9o+n ~ l /L?. By Chow's theorem [9], 
Prym(C/Co) is thus a subabelian variety of Jac (C) . More precisely, we have 

Prym(C/C ) = (KerN v )°, 

= Ker(l Jac{c) + cr)°, 

= Im { l Jac{c) - o-) C Jac(C). 

Equivalently, the involution a induces an involution 

a : Jac(C) — ► Jac(C), class of D i — > class of aD, 

and up to some points of order two, Jac(C) splits into an even part Jac (Co) 
and an odd part Prym(C/Co) : 



Jac(C) = Prym{C/C ) Jac (C ) , 
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with 



dimJac(Co) = go, 

dim Jac() = g = 2go + n — 1, 
dim Prym(C /Co) = 5 - 50 = 50 + n - 1. 

Observe that A (2) (resp. T(3)) is the period matrix of Jac(Co) (resp. 
Prym(C/C )). Write 

r = (i/, v) , 



with 



/ 2/ a 



2 f„ L w i • • • f„ w i \ 



/ 2 4^1 



2 /b w l 

J0 90 



/I ^1 \ 



V 



and let us call 

/ 2 Ja^l X 



ei 



J6 so+ 1 

2 Ib 90 Ugo+n-l J bg(j+1 Ug +n-l ... J 6go+n _ 1 W po+ „-i y 
/ 2 L Wi \ 



\ 2 Jai W 90+n-l / 



e 90+l 



V 2 /a go / 



\ -^so+i / 

Then, in the new basis (Ai, . . . , A ff0+n _i) where 



1 pour 1 < j < g , 

2 pour 50 + I <j < go + n - 1, 



the period matrix T takes the canonical form 

(A 5 , Z), 

with 



and 



A s = diag(5i,. ..,S n ) 
Z = AgU^V, 
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symmetric and ImZ > 0. Then 

r * = (^SO+n-lAjS^o+n-lAj^Aj 1 ) , 

and 

r* = (u*,v*), 

= (A B D E), 

is the period matrix of the dual abelian variety Prym* (C / Co) ■ The above 
discussion is summed up in the following statement : 

Theorem 1 Let 

<p:C^Co, 

be a double covering where C and Cq are nonsingular algebraic curves with 
jacobians Jac(C) and Jac(Co). Let 

a : C — ► C, 

be the involution exchanging sheets of C over Co = C/a. This involution 
extends by linearity to a map (which will again be denoted by a) 

a : Jac{C) — ► Jac(C), 

and up some points of order two, Jac{C) splits into an even part ,i.e., Jac{C) 
and an odd part (called Prym variety) denoted Prym(C /Cq) and defined by 

Prym(C/C ) = (h° (C, Sift ") * /H, (C, Z)" , 

where is the sheaf of holomorphic 1- forms on C and — denote the -1 
eigenspace for a vector space on which j acts. To be precise, we have 



with 



and 



Jac(C) = Prym(C/C ) Jac (Co) , 

dim Jac (Co) = genus go of Co, 
dim Jac (C) = genus g of C = 2go + n — 1, 



dim Prym(C/C ) = g - g = go + n + 1, 
with 2n branch points. Moreover, if 



ft 



A B C D E F 
G H I J K L 
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is the period matrix of Jac{C) where A, . . . ,L denote the matrices (1) , then 
the period matrices of Jac(Co) , Prym {C/Cq) and the dual Prym* (C/Co) of 
Prym(C/Co), are respectively 

A = (G H), 
r = (2A B 2D E), 

and 

r* = (A B D E). 

3 Algebraic complete integrability 

3.1 A survey on abelian varieties and algebraic integrability 

We give some basic facts about integrable hamiltonian systems and some 
results about abelian surfaces which will be used, as well as the basic tech- 
niques to study two-dimensional algebraic completely integrable systems. 
Let T = C/A be a n— dimensional abelian variety where A is the lattice 
generated by the 2n columns Ai, . . . , \2n of the n x 2n period matrix O and 
let D be a divisor on T. Define 

C{V) = {/ meromorphic on T : (/) > -V}, 

i.e., for V = kjT>j a function / € £(D) has at worst a kj— fold pole along 
Vj. The divisor T> is called ample when a basis (/o, . . . , /at) of C(kT>) embeds 
T smoothly into for some k, via the map 

T^F N , p— [l:/i(p):... : f N (p)], 

then kT> is called very ample. It is known that every positive divisor V on 
an irreducible abelian variety is ample and thus some multiple of V embeds 
M into F N . By a theorem of Lefschetz, any k > 3 will work. Moreover, there 
exists a complex basis of C n such that the lattice expressed in that basis is 
generated by the columns of the n x 2n period matrix 

(Si | \ 

| z . 

V o s n | / 

with Z T = Z,lmZ > 0,Sj £ N* and Sj\Sj + \. The integers Sj which provide 
the so-called polarization of the abelian variety M are then related to the 
divisor as follows : 

dimC(V) = 5 1 ...5 n . (4) 
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In the case of a 2— dimensional abelian varieties (surfaces), even more can be 
stated : the geometric genus g of a positive divisor V (containing possibly 
one or several curves) on a surface T is given by the adjunction formula 

, x K T .V + V.V 

giP) = + 1, (5) 

where Kt is the canonical divisor on T, i.e., the zero-locus of a holomorphic 
2— form, V.V denote the number of intersection points of V with a+V (where 
a + V is a small translation by a of V on T), where as the Riemann-Roch 
theorem for line bundles on a surface tells you that 

X {V)=p a {T) + l+ l -(V.V-VK M ), (6) 

where p a (T) is the arithmetic genus of T and the Euler characteristic 

of V. To study abelian surfaces using Riemann surfaces on these surfaces, 
we recall that 

X (V) = dim H (T,O T (V)) - dim H\T,O t (D)), 

= dim£(V) - dim if 1 (T, Q 2 (V ® K%)), (Kodaira-Serre duality), 
= dim C(V), (Kodaira vanishing theorem), (7) 

whenever T>®K^ defines a positive line bundle. However for abelian surfaces, 
Kt is trivial and p a (T) = —1; therefore combining relations (4), (5), (6) and 
(7), we obtain 

X (V) = dim £{V), 
V.V 



2 ' 

= g(p)-l, 

= 5^2. (8) 

Let M be a 2n-dimensional differentiable manifold and u a closed non- 
degenerate differential 2-form. The pair (M, to) is called a symplectic man- 
ifold. Let H : M — ► R be a smooth function. A hamiltonian system on 
(M, u) with hamiltonian H can be written in the form 

dH dH dH dH 

ill = -K— ,-,qn = -7T- , Pi = -"5— ,-,Pn = , (9) 

dpi dp n dqi dq n 

where (qi, q n ,Pi, ■■■,Pn) are coordinates in M. Thus the hamiltonian vec- 
tor field Xh is defined by 

^ f dH d dH d 



k=l 



dpk dq k dq k dp k 
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If F is a smooth function on the manifold M, the Poisson bracket {F, H} of 
F and H is defined by 



X H F 



dH dF dH dF 



k=i 
{F,H} 



dp k dq k dq k dp k 



(10) 



A function F is an invariant (first integral) of the hamiltonian system (9) if 
and only if the Lie derivative of F with respect Xh is identically zero. The 
functions F and H are said to be in involution or to commute, if {F, H} = 0. 
Note that equations (9) and (10) can be written in more compact form 

dH 

± = J ~dx' x = {<lli--i ( ln,Pl,~;Pn) T , 



{F,H} 



ldF_ jdH\ 

\dx' dx I ' 



k,l=l 



dF dH 

dx k dxi ' 



with J = 
matrix anc 



O 



I 

o 



, a skew-symmetric matrix where I is the n x n unit 



O the n x n zero matrix. A hamiltonian system is completely inte- 
grable in the sense of Liouville if there exist n invariants Hi = H, H2, ■ ■ ■ , H n 
in involution (i.e., such that the associated Poisson bracket {H k ,Hi} all van- 
ish) with linearly independent gradients (i.e., dH\ A ... A dH n 7^ 0). For 
generic c = (ci, c n ) the level set 

M c = {Hi = d,...,H n = c n }, 

will be an n-manifold, and since 

x H k Hi = {H k , Hi} = 0, 

the integral curves of each Xff fc will lie in M c and the vector fields Xu k span 
the tangent space of M c . By a theorem of Arnold [4,15], if M c is compact 
and connected, it is diffeomorphic to an n-dimensional real torus and each 
vector field will define a linear flow there. To be precise, in some open 
neighbourhood of the torus one can introduce regular symplectic coordinates 



si, . 



in which uj takes the canonical form 



to = ds k A dip k . 



k=i 

Here the functions s k (called action- variables) give coordinates in the direc- 
tion transverse to the torus and can be expressed functionally in terms of 
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the first integrals H k . The functions <p k (called angle- variables) give standard 
angular coordinates on the torus, and every vector field Xjj k can be written 
in the form 

= h k (si, . . . ,s n ) , 

that is, its integral trajectories define a conditionally-periodic motion on the 
torus. Consequently, in a neighbourhood of the torus the hamiltonian vector 
field Xn k take the following form 

4 = 0, (p k = h k (si,...,s n ), 

and can be solved by quadratures. 

Consider now hamiltonian problems of the form 

X H :± = J— = /(*), x G R m , (11) 

where H is the hamiltonian and J = J(x) is a skew-symmetric matrix with 
polynomial entries in x, for which the corresponding Poisson bracket 

satisfies the Jacobi identities. The system (11) with polynomial right hand 
side will be called algebraically completely integrable (a.c.i.) in the sense of 
Adler-van Moerbeke [1] when : 

(i) The system possesses n+k independent polynomial invariants Hi, H n+k 
of which k lead to zero vector fields 

jdHn+i l<i<k, 
ox 

the n remaining ones are in involution (i.e., {Hi,Hj} = 0) and m = In + k. 

n+k 

For most values of q G R, the invariant varieties f] {x £ R m : Hi = q} are 

i=i 

assumed compact and connected. Then, according to the Arnold-Liouville 
theorem, there exists a diffeomorphism 

n+k 

f]{x£ R m : Hi = ci} -> R n /Lattice, 
i=i 

and the solutions of the system (11) are straight lines motions on these tori. 

(ii) The invariant varieties, thought of as affine varieties 

n+k 

M c = f]{Hi = a,x eC ra }, 
i=i 

in C m can be completed into complex algebraic tori, i.e., 

M c U V = C 1 / Lattice, 

where C n / Lattice is a complex algebraic torus (i.e., abelian variety) and V 
a divisor. 
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Remark 3.1 Algebraic means that the torus can be defined as an intersec- 
tion 

f]{P l (X ,...,X N )=0}, 

i 

involving a large number of homogeneous polynomials Pi. In the natural 
coordinates (ti,...,t n ) of C n /Lattice coming from C n , the functions Xi = 
Xi(ti, t n ) are meromorphic and (11) defines straight line motion on C n / Lattice. 
Condition (i) means, in particular, there is an algebraic map 

(xi(t),...,X m {t)) I > (fj,i(t),...,fl n (t)), 

making the following sums linear in t : 




= djt , 1 < j < n, dj G C, 



where wi,...,w„ denote holomorphic differentials on some algebraic curves. 

Adler and van Moerbeke [1] have shown that the existence of a coherent 
set of Laurent solutions : 

oo 

Xi = ^Tx ( : i) t j - ki , ki G Z, some h > 0, (12) 

3=0 

depending on dim (phase space) — 1 = m — 1 free parameters is necessary 
and sufficient for a hamiltonian system with the right number of constants 
of motion to be a.c.i. So, if the hamiltonian flow (11) is a.c.i., it means 
that the variables Xi are meromorphic on the torus C n / Lattice and by com- 
pactness they must blow up along a codimension one subvariety (a divisor) 
V C C™ / 'Lattice. By the a.c.i. definition, the flow (11) is a straight line 
motion in C n / Lattice and thus it must hit the divisor V in at least one 
place. Moreover through every point of V, there is a straight line motion 
and therefore a Laurent expansion around that point of intersection. Hence 
the differential equations must admit Laurent expansions which depend on 
the n — 1 parameters defining V and the n + k constants Cj defining the torus 
C n /Lattice , the total count is therefore m — 1 = dim (phase space) — 1 
parameters. 

Assume now hamiltonian flows to be (weight)-homogeneous with a weight 
going with each variable Xj, i.e., 

f i (a ul x 1 ,...,a Um x m ) =a Vi+1 f i (x 1 ,...,x m ), Vq G C. 

Observe that then the constants of the motion H can be chosen to be 
(weight)-homogeneous : 

H (a Ul xi, ...,a Um x m ) = a k H (xi, ...,x m ) , k G Z. 
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If the flow is algebraically completely integrable, the differential equations 
(11) must admits Laurent series solutions (12) depending on m — 1 free 
parameters. We must have k{ = vi and coefficients in the series must satisfy 
at the O t/l step non-linear equations, 

fi {xf\...,x^)+ gi xf ] =0, 1 < i < m, (13) 
and at the k th step, linear systems of equations : 

(L-kI)z^ = { , . 1 .° % k = \ k 1)f h 1 (14) 

[ some polynomial m x v x y 1 tor k > 1, 

where 

L = Jacobian map of (13) , 
df 

If m — 1 free parameters are to appear in the Laurent series, they must ei- 
ther come from the non- linear equations (13) or from the eigenvalue problem 
(14) , i.e., L must have at least m — 1 integer eigenvalues. These are much 
less conditions than expected, because of the fact that the homogeneity k 
of the constant H must be an eigenvalue of L Moreover the formal series 
solutions are convergent as a consequence of the majorant method [1] . So, 
the question is how does one prove directly that the system is effectively 
a.c.i. with abelian space coordinates? The idea of the direct proof used 
by Adler and van Moerbeke is closely related to the geometric spirit of the 
(real) Arnold-Liouville theorem discuted above. Namely, a compact complex 
n-dimensional variety on which there exist n holomorphic commuting vector 
fields which are independent at every point is analytically isomorphic to a 
n-dimensional complex torus C n / Lattice and the complex flows generated 
by the vector fields are straight lines on this complex torus. Now, the affine 
variety M c is not compact and the main problem will be to complete M c into 
a non singular compact complex algebraic variety M c = M c U V in such a 
way that the vector fields Xh x , Xn m extend holomorphically along the di- 
visor V and remain independent there. If this is possible, M c is an algebraic 
complex torus, i.e., an abelian variety and the coordinates Xi restricted to 
M c are abelian functions. A naive guess would be to take the natural com- 
pactification M c of M c in P m (C). Indeed, this can never work for a general 
reason: an abelian variety M c of dimension bigger or equal than two is never 
a complete intersection, that is it can never be descriped in some projective 
space P m (C) by m-dim M c global polynomial homogeneous equations. In 
other words, if M c is to be the affine part of an abelian variety, M c must 
have a singularity somewhere along the locus at infinity M c n {Xq = 0}. In 
fact, Adler and van Moerbeke [1] showed that the existence of meromorphic 
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solutions to the differential equations (2) depending on m — 1 free param- 
eters can be used to manufacture the tori, without ever going through the 
delicate procedure of blowing up and down. Information about the tori can 
then be gathered from the divisor. An exposition of such methods and their 
applications can be found in [1,2]. 



3.2 The Henon-Heiles system 

The Henon-Heiles system 

dH dH dH dH 

dpi dp 2 oqi dq 2 

with 

H = Hi = - (p\ + p\ + aq\ + bq 2 2 ) + q\q 2 + 6q 3 , 

has another constant of motion 

H 2 = qf + kq\q\ - 4pi (piq 2 - p 2 qi) + ^aqfq 2 + (4a - b) (pf + aqf) , 

where a, b, are constant parameters and qi,q 2 ,Pi,p 2 are canonical coordi- 
nates and momenta, respectively. First studied as a mathematical model 
to describe the chaotic motion of a test star in an axisymmetric galactic 
mean gravitational field this system is widely explored in other branches of 
physics. It well-known from applications in stellar dynamics, statistical me- 
chanics and quantum mechanics. It provides a model for the oscillations of 
atoms in a three-atomic molecule. The system (15) possesses Laurent series 
solutions depending on 3 free parameters a, (3, 7, namely 

.1 = J+(^ + f-^)t + ^ + q^ + qf¥ + q r¥ + ..., 
1 t a 2 B t fa 4 a 2 A a 2 B B 2 \ ^ ( a(3 3 4 



92 = -^ + 12-12 + U + ^""60-" 240r + "3"^ + 7t + '" 
where Pi = q±, p 2 = q 2 and 

( 4 ) aAB a 5 lla 3 B lla 3 A aB 2 aA 2 
Ql ~ 24 ~ 72 + 720 120 "720~ 8~' 

(5) _ Pa 2 0B A0 
qi 12 60 10 ' 

(6 ) 07 a 7 a 5 A a b B a 3 B 2 a 3 A 2 a 3 AB aAB 2 
Ql 9~ ~ 15552 ~ 2160 + 12960 + 25920 + 1440 ~ 4320 + 1440 

aB 3 aA 2 B aA 3 



+ 



19440 288 144 
Let V be the pole solutions restricted to the surface 
2 

M c = n { x = (?l'92,Pl,P2) G C 4 ,Hi (x) = Ci] 

1=1 
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to be precise T> is the closure of the continuous components of the set of 
Laurent series solutions x (t) such that Hi(x(t)) = Ci, 1 < i < 2, i.e., 
V = t — coefficient of M c . Thus we find an algebraic curve defined by 

V:[3 2 = P 8 (a), (16) 

where 

„ , x 7 o 1 / „ 13 \ fi 1 / 671 2 17 ,n 943 „ „ \ , 

P s ( a ) = a 8 5A B a 6 B 2 + —A 2 BA a 4 

8 v ' 15552 432 V 18 / 36 V 15120 7 1260 1 



1 ( . ,s 1 «s 13 ,, n 2 „ o 10 \ n 1 

AA 3 B 3 A 2 B + -AB 2 c x a 2 + — c 2 . 

36 V 2520 6 9 7 V 36 

The curve T> determined by an eight-order equation is smooth, hyperelliptic 
and its genus is 3. Moreover, the map 

a:V^V, (J3,a)^>(p,-a), (17) 

is an involution on V and the quotient 8 = V/a is an elliptic curve defined 
by 

£ : P 2 = P 4 (C), (18) 

where P4 (£) is the degree 4 polynomial in C = a 2 obtained from (16) . The 
hyperelliptic curve T> is thus a 2-sheeted ramified covering of the elliptic 
curve £ (18) , 

p:V^£, (fta) — (/J,C), (19) 

ramified at the four points covering C = and 00. Following the methods 
in Adler-van Moerbeke [1,2], the affme surface M c completes into an abelian 
surface M c , by adjoining the divisor T>. The latter defines on M c a polariza- 
tion (1,2). The divisor 2T> is very ample and the functions 

1, yi, Vi, 2/2, zi, x 2 + y 2 y 2 , x 2 yi-2xiy 2 , x 1 x 2 + 2Ay 1 y 2 + 2y 1 y 2 , 

embed M c smoothly into CP 7 with polarization (2, 4), according to (8). Then 
the system (15) is algebraic complete integrable and the corresponding flow 
evolues on an abelian surface M c = C 2 /lattice, where the lattice is generated 

uxi. ■ 1 f 2 a c \ t f a c \ n 

by the period matrix I 1 , 1m I I > U. 

Theorem 2 T/ie abelian surface M c which completes the affine surface M c 
is the dual Prym variety Prym* (T>/£) of the genus 3 hyperelliptic curve V 
(16) for the involution a interchanging the sheets of the double covering p 
(19) and the problem linearizes on this variety. 
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Proof. Let (ai, a 2 , a 3 , 61, 62 , ^3) be a canonical homology basis of V such that 

0-(ai) = a 3 , <t(6i) = 6 3 , a (a 2 ) = -a 2 , cr(& 2 ) = -& 2 , 

for the involution u (17). As a basis of holomorphic differentials wo , wi , ^2 
on the curve V (16) we take the differentials 

a 2 da da ada 



ui = —5—, w 2 = — , lo 3 



and obviously 

£7*(o>i) = -wi, ct*(cj 2 ) = -W2, o-*(u s ) = 0J 3 . 
We consider the period matrix Q of Jac(T>) 

I 4 ^1 4 4 w i 4 w i 4 w i 4 w i 

n = 4 ^2 4 ^2 4 ^ 2 4 ^2 4 ^2 4 ^2 

V 4 ^3 4 ^3 4 ^3 4 ^3 4 ^3 4 ^3 

By theorem 1, 

/ 4 Wi 4 Wi -4 Wi 4 Wi 4 Wi -4 Wi 
n = 4 ^ 2 4 ^ 2 " 4 ^ 2 4 ^ 2 4 ^ 2 - 4 ^ 2 
V 4 ^ 3 4 ^ 3 4 W3 4 ^ 3 

and therefore the period matrices of Jac{£) (i.e., £), Prym(V / £) and Prym* (V / '£) 
are respectively 

A = ( / ^3 / W 3 ), 

r = f 2 ^1 Wl ^« Wl 2 ^ Wl h Ul 

V 2 4^ 2 4 W2 2 4^ 2 4^ 2 

and 



Let 



T * = ( 4 w i 4 Wi 4 Ui 4 Wi 

v 4 U2 4 W2 4 ^ 2 4 W2 

Ln = (^2m,i uj 2 + / w 2 : m^n; € Z}, 

i=i V w 3 / A< \ "3 / 
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be the period lattice associated to f2. Let us denote also by La, the period 
lattice associated A. According to theorem 1, we get the following diagram 



I 

£ V 

— ► kerA^ — ► Prym{V / '£) © £ = Jac{V) ^ £ — > 
\r _ | 

M c = M c U 2V ~ C 2 /lattice 

I 


The polarization map 

r : Prym(V/£) — > M c = Prym*{V/£), 

has kernel (<f*£) — 1>2 x ^2 and the induced polarization on Prym(D / '£) is 
of type (1,2). Let 



be the uniformizing map where dt\,dt2 are two differentials on M c corre- 
sponding to the flows generated respectively by H±,H2 such that : 



dti\x> = wi, dt2\v = 

and 

L A = {J2n k ( j"** 1 J :ii fe £Z}, 

fc=i V Jffc 2 / 

is the lattice associated to the period matrix 

j^i ;„ 4 *i l 4 ^i 
j ui dt 2 j U2 dt 2 j Ki dt 2 j U4 dt 2 

where (v\, z^, ^3, ^4) is a basis of Hi(M c ,Z). By the Lefschetz theorem on 
hyperplane section [9], the map 

#i(2>,Z) — £Ti(M c ,Z), 

induced by the inclusion P ^ M c is surjective and consequently we can find 
4 cycles 1/1, ^2, ^3, ^4 on the curve D such that 

A = 1 L, "1 L 2 ^1 L 4 L 4 

1^2 f V2 ^2 J U3 U 2 f U4 UJ 2 
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and 

L A = {)^n k I I :n k eZ}. 

k =i V J "* 2 / 

The cycles ^ljj^ , ^ 3 , ^4 in ^ which we look for are 01,61,02,62 and they 
generate Hi(M c ,Z) such that 




A 



/ 01 wi / 6l wi J^i / fe2 ^i 

/ ai W2 / 6l W2 J a2 ^2 



is a Riemann matrix. We show that A = T* ,i.e^, the period matrix of 
Prym*(T>/£) dual of Prym(D/£). Consequently M c and Prym*(D/£) are 
two abelian varieties analytically isomorphic to the same complex torus 
C 2 /L\. By Chow's theorem, M c and Prym*(V/£) are then algebraically 
isomorphic. 

3.3 The Kowalewski rigid body motion 

The motion for the Kowalewski's top is governed by the equations 

m = m A Am + 7 A I, (20) 
7 = 7 A Am, 

where m,7 and I denote respectively the angular momentum, the direc- 
tional cosine of the z-axis (fixed in space), the center of gravity which 
after some rescaling and normalization may be taken as / = (1,0,0) and 
Am = (mi/2, m-z/2, m 3 /2) . The system (20) can be written 

mi = m 2 m 3 , 

m 2 = -mim.3 + 273, 

m 3 = -2 72 , (21) 
7! = 2m 3 72 -m 2 73, 

72 = mi7 3 -2m 3 7i, 

73 = m 27i - "H72, 

with constants of motion 

Hi = - {m\ + m|) + m\ + 271 = ci, 

iJ 2 = mi7i + m 2 7 2 + m 3 7 3 = c 2 , (22) 
#3 = 7i 2 + 7 2 2 + 7l = c 3 = 1, 

H i = [ ( mi ^ Z7W2 ) - (71 + ^2) I [ ( mi ) - (7i - *7!2) I = c 4- 
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The system (21) admits two distinct families of Laurent series solutions 

mim _J f+*(«?-2)a 2 + o(t), 
1 W " I t- t (a?-2)a 2 + o(t), 



m 2 (i) = 
m 3 (t) = 
71 (*) 
72(0 
73 (t) 



m _ a j a2 + o (t) , 
=fi - af a 2 + o (i) 

| + aia 2 + o (t) , 
=f + aia 2 + o (t) , 

2^ + ° (*) , 

2F + ° (*) . 

# + ° (t) , 

f +o(t), 



which depend on 5 free parameters ai, a$. By substituting these series in 
the constants of the motion Hi (22), one eliminates three parameters linearly, 
leading to algebraic relation between the two remaining parameters, which 
is nothing but the equation of the divisor V along which the rrii,^ blow up. 
Since the system (21) admits two families of Laurent solutions, then V is a 
set of two isomorphic curves of genus 3, V = V\ + T>-\ : 

V £ : P (ai, a 2 ) = (a? - l) ((a? - l) a\ - P (a 2 )) + c 4 = 0, (23) 

where 

P (a 2 ) = cia 2 - 2ec 2 a 2 - 1, 

and e = ±1. Each of the curve P e is a 2 — 1 ramified cover (ai,a 2 ,/3) of 
elliptic curves T>® : 

V° £ : /3 2 = P 2 (a 2 ) - 4c 4 a|, (24) 

ramified at the 4 points «i = covering the 4 roots of P (a 2 ) = 0. It 
was shown by the author [14] that each divisor V £ is ample and defines a 
polarization (1, 2) , whereas the divisor D, of geometric genus 9, is very ample 
and defines a polarization (2,4) , accordind to (8). The affine surface 

4 

M c = P| {Hi = c t } C C 6 , 
1=1 

defined by putting the four invariants (22) of the Kowalewski flow (21) equal 
to generic constants, is the affine part of an abelian surface M c with 

M c \ M c = T> = one genus 9 curve consisting of two genus 3 
curves T> £ (23) intersecting in 4 points. Each 
T> £ is a double cover of an elliptic curve P° (24) 
ramified at 4 points. 



24 



Moreover, the Hamiltonian flows generated by the vector fields Xjjx and Xff 4 
are straight lines on M c . The 8 functions 

1, fi = m 1 , f2 = m 2 , /3 = m 3 , / 4 = 73, h = fi+fh 

h = 4 /i/4 - /3/5, /7 = (/271 - /172) /3 + 2/472, 

form a basis of the vector space of meromorphic functions on M c with at 
worst a simple pole along V Moreover, the map 

M c ~ C 2 /Lattice - CP 7 , (ti,* 2 ) h-» [(1, / x (t l5 1 2 ) , .... / 7 (ti, t 2 ))] , 

is an embedding of M c into CP . Following the method (subsection 3.1), we 
obtain the following theorem : 

Theorem 3 The tori M c can be identified as 

M c = Prym*(V £ /V° £ ), 
i.e., dual of Prym(T> £ /'D £ ) and the problem linearizes on this Prym variety. 

3.4 Kirchhoff's equations of motion of a solid in an ideal fluid 

The Kirchhoff's equations of motion of a solid in an ideal fluid have the form 



Pi 


dH 

oh 


dH 

dl 2 


I 






dH 
oh 


dH 

oh 


t 




P3 


dH 

dh 


dH 

oh 


) 




h 


dH 

dp 3 


dH 

dp 2 


, dH 

dl 3 


. dH 
~ h dh 


12 


dH 

dpi 


dH 

dp 3 


, dH 
dh 


. dH 
~ h dh 


h 


dH 

= P 1 ^ 
dp2 


dH 
dpi 


, dH 

dh 


. dH 
dh 



(25) 



where (pi,P2,P3) is the velocity of a point fixed relatively to the solid, 
(^l) h, h) the angular velocity of the body expressed with regard to a frame 
of reference also fixed relatively to the solid and H is the hamiltonian. 
These equations can be regarded as the equations of the geodesies of the 
right-invariant metric on the group E (3) = SO (3) x R 3 of motions of 3- 
dimensional euclidean space M 3 , generated by rotations and translations. 
Hence the motion has the trivial coadjoint orbit invariants (p,p) and (p,l). 
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As it turns out, this is a special case of a more general system of equations 
written as 

dH dH 

x = x A — + yA —, 
ox dy 

dH dH 

y = y a — + x a — , 

ox oy 

where x = (xi,X2,xs) £ M 3 et y = (2/1,1/2,2/3) G I^ 3 - The first set can be 
obtained from the second by putting (x,y) = (l,p/e) and letting e — > 0. 
The latter set of equations is the geodesic flow on 50(4) for a left invariant 
metric defined by the quadratic form H. In Clebsch's case, equations (25) 
have the four invariants : 

Hi = H = - (a lP j + a 2 p 2 2 + a 3 p 2 3 + b x l\ + b 2 l\ + 6 3 /|) , 

H 2 = pl+pl+pl (26) 

^3 = Pih+P2h+P3k, 



with 



Ha = - {b 1 pl + b 2 pl + hpl + q{1\ + + 
a.2 ~ a 3 + Q3 ~ Qi + Qi ~ Q-2 _ Q 



h b 2 b 3 

and the constant g satisfies the conditions 

_ h (6 2 ~ 63) _ h (b 3 - h) _ 63 (61 - 6 2 ) 
a 2 - 03 «3 - «i ai - a,2 

The system (25) can be written in the form (11) with m = 6; to be precise 

dH 

x = f(x) = J— , x = (pi,P2,P3,h,l2,h) J , (27) 

where 



P = 



L 





P 




-( 


P L 


)■ 





~P3 


P2 


P3 





-P\ 




Pi 








~h 




^3 





-il 


-«2 


Pi 






Consider points at infinity which are limit points of trajectories of the flow. 
In fact, there is a Laurent decomposition of such asymptotic solutions, 

x(t)=t- 1 ^+x^t + x^t 2 + ...), (28) 
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which depend on dim(phase space) — 1 = 5 free parameters. Putting (28) into 
(27), solving inductively for the x^ k \ one finds at the th step a non-linear 
equation, 

x (o) + /(x (o) ) = 0) 
and at the k th step, a linear system of equations, 

(L - kI) X W = I for fc = 1 

\ quadratic polynomial in x^ l \ for k > 1, 

where L denotes the jacobian map of the non-linear equation above. One 
parameter appear at the th step, i.e., in the resolution of the non-linear 
equation and the 4 remaining ones at the k th step, k = 1, ...,4. Taking into 
account only solutions trajectories lying on the invariant surface 

4 

M c = P| {Hi (x) = a} C C 6 , 
1=1 

we obtain one-parameter families which are parameterized by a curve. To be 
precise we search for the set of Laurent solutions which remain confined to 
a fixed afflne invariant surface, related to specific values of c±, 02,03,04, i.e., 

4 

V = (~){t° - coefficient of Hi (x (t)) = a} , 
i=\ 

= an algebraic curve defined by 

9 2 + ci/3 2 7 2 + c 2 aV + c 3 a 2 /3 2 + c 4 a^ = 0, (29) 

where 6 is an arbitrary parameter and where 

(0) a (°) (°) 

parameterizes the elliptic curve 

£ : P 2 = d\a 2 - 1, 7 2 = a%a 2 + 1, (30) 
with di,d,2 such that: 

d\ + a% + 1 = 0. 

The curve V is a 2-sheeted ramified covering of the elliptic curve £. The 
branch points are defined by the 16 zeroes of ci/3 2 7 2 +C2a 2 7 2 +C3a 2 /3 2 +C4a/37 
on £. The curve T> is unramified at infinity and by Hurwitz's formula, the 
genus of V is 9. Upon putting £ = a 2 , the curve V can also be seen as a 
4— sheeted unramified covering of the following curve of genus 3 : 

C : (6 2 + Cl /3 2 7 2 + (c 27 2 + c 3 (3 2 ) C) 2 - c 2 (f3 2 7 2 = 0. 
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Moreover, the map 

r: C — C, (0,0 — (-0,0, 

is an involution on C and the quotient Co = C/r is an elliptic curve defined 
by 

Co : v 2 = 4C (dl4( 2 + (d\ -4)(-l). 
The curve C is a double ramified covering of Co, 

c — Co, (M.C)-- fa.C), 

fl 2 = - Cl /3 2 7 2 - (c 27 2 + c 3 /3 2 ) C + V 

v 2 = cjc (44c 2 + (4 - 4) C - 1) ■ 

Let (a±, a 2 , CI3, 61, 62, 63) be a canonical homology basis of C such that 



C 



r(ai) = a 3 , r(6i) = 6 3 , r(o 2 ) = -a 2 , r(6 2 ) = -6 2 , 
for the involution r. Using the Poincare residu map [9], we show that 

d( (d( d( 

77 0ry 6/77 

form a basis of holomorphic differentials on C and 

t*(uj )=uj , t* (u k ) = -u k , (A; = 1,2). 

Haine [10] shows that the flow evolues on an abelian surface M c C CP 7 of 
period matrix ( ° ^ ] , Im ( a ? ) > and also identified M c as 



Prym variety Prym{C / Cq) . 

Theorem 4 27te abelian surface M c can be identified as Prym{C / Cq) . More 
precisely 

4 

P| {x e C 6 , Hi (x) = a} = Prym(C/C )\V, 
i=i 

where V is a genus 9 curve (29), which is a ramified cover of an elliptic 
curve £ (30) with 16 branch points. 
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